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abstract 

We have examined the distribution of HI and H 2 gases in four face-on galaxies by using 
the observed data of CO and HI line emissions from the literatures. We demonstrate that 
the gas phase transition of HI and H 2 occurs suddenly within a narrow range of radius, 
which we call the molecular front. We have tried to explain such phase transition in 
galactic scale with a help of the phase transition theory proposed by Elmegreen. The 
crucial parameters for determinating the molecular fraction f mo \ are interstellar pressure 
P , UV radiation held U, and metallicity Z, and we have constructed a model galaxy in 
which P , U and Z obey an exponential function of the galacto-centric radius. The model 
shows that the molecular front must be a fundamental feature of galaxies which has an 
exponential disk, and that the metallicity gradient is most crucial for the formation of 
the molecular front. We have also tried to reproduce the observed molecular fraction f mo \ 
by giving the set of (P, U, Z) observationally, and show that the model can describe the 
variation of the molecular fraction f mo \ in galaxies quite well. We discuss the implication 
of the molecular front for the chemical evolution of galaxies. 


keywords: galaxies - atomic gas - molecular gas - gas phase transition - evolution of 
galaxies 
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1 Introduction 


The HI and H 2 gases are the main components of the interstellar matter (ISM) and 
both are tightly related to the star forming activity in galaxies. Therefore, investigating 
these neutral components of the ISM is one of the most powerful ways to understand 
the structure and evolution of galaxies. The HI gas can be observed directly in the HI- 
21cm line emission and H 2 can be investigated by observing the CO-2.6mm line, while 
it is difficult to observe H 2 molecules directly. So far we have observed many galaxies 
in the CO-2.6mm line emission using the Nobeyama 45m telescope (e.g. Sofue & Nakai 
1993; 1994; Sofue 1994). Many other researchers have observed galaxies in the HI-21cm 
line using VLA and other facilities (e.g. Rupen 1991; Irwin & Seaquist 1991). However, 
in spite of a large number of observations of galaxies in CO and HI, there has been 
yet no direct comparison of the distribution of HI and H 2 with sufficiently high spatial 
resolutions. 

We have made the first detailed comparison of the distributions of HI and H 2 (calculated 
from CO) in several edge-on galaxies (Sofue, Honma & Arimoto 1994, hereafter Paper 
I) based on derivations of the density distributions of HI and H 2 . We have obtained the 
radial variation of the molecular fraction (/ mo i) which is defined by the ratio of molecular 
gas density to total hydrogen gas density as, 


/mol 


2 x n(H 2 ) 


n(H) + 2 x n(H 2 ) ’ 


( 1 ) 


By making such quantitative comparison, we have demonstrated 1) that the central region 
within a few kpc is totally molecular, 2) that the outer region is totally atomic, and 3) 
that the HI-H 2 gas phase transition occurs rather sharply at a certain radius, which we 
shall call the molecular front. The disk of a galaxy is divided into two radial zones at the 
molecular front; the inner H 2 disk and the outer HI disk. 

On the other hand, a theoretical study on the phase transition between HI and H 2 has 
been made by Hollenbach et al.(1971) in rather small scale corresponding to individual 
sizes of clouds and has recently been developed by Elmegreen (1993) in galactic scale. 
However, there has been no study yet that argues how / mo i varies with the distance from 
the galactic center or whether a feature like the molecular front can be theoretically 
reproduced. In this paper we construct a gas phase transition model for disk galaxies 
and investigate the variation of / mo i with the galacto-centric radius. We show that the 
observed features mentioned above, including the molecular front, can be reproduced by 
the model. 


2 The Molecular Front in Face-on Galaxies 

In Paper I, we obtained the distributions of HI and H 2 gases in four edge-on galaxies by 
deconvolving position-velocity diagrams under the assumptions of a flat rotation curve and 
an axisymmetric distribution. In order to avoid the uncertainties which might come from 
these assumptions, we now investigate the distributions of HI and H 2 gases in nearby 
face-on disk galaxies and see whether the feature like the molecular front claimed in 
Paper I really appears in face-on galaxies. From many face-on galaxies so have been 
observed in HI and CO, we select a sample of four galaxies which are near enough to get 
sufficiently high resolution and are investigated well to obtain other properties such as 
the surface brightness and the metallicity gradient. Our sample consists of four galaxies; 
M51, M101, NGC 6946 and IC342, which are typical late-type spirals classified as Sc, and 
many fundamental observational data are available including HI, CO, optical luminosity 
profile and oxygen abundance of HII regions. The parameters for these galaxies and the 
references are summarized in Table 1. 

Fig. l.a - l.d 
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We show the surface density distributions of Hf and H 2 gases in these four galaxies in 
Fig.l.a - l.d, where the amount of H 2 is calculated from CO intensity using a constant 
conversion factor of 3xl0 20 (H 2 (K km/s) -1 ). We conventionally assume that the con¬ 
version factor is universal and is the same as the Galactic value. The gas distribution for 
M101 has been azimuthally averaged, while those for other galaxies have been obtained 
by averaging the distributions along the major axis in both sides. All galaxies in Fig.l 
show conspicuous concentration of molecular gas at the galactic center with the surface 
mass density of more than 100 M©/pc 2 . The density of molecular gas decreases steeply 
with the radius, showing local enhancement in the arms in the cases of M51 and IC 342. 
On the other hand, the distribution of HI gas is rather flat, almost constant through the 
galactic disk. The enhancement of HI density in the arm is not so clear as seen in the H 2 . 
This implies that the gas might be converted into H 2 once the gas density becomes larger 
than the critical density which is around 5 M 0 /pc 2 . 

Fig, 2.a 2,b 

In Fig.2.a we show the molecular fraction f mo \ which is defined by equation (1). We can 
see that f mo i at the galactic center is almost unity in any of these galaxies and that f mo \ 
decreases drastically at a certain radius, except for M51. For M51, the molecular fraction 
is still high at the maximum radius of our data, suggesting that the phase transition must 
occur at a larger radius. The characteristic behavior of / mo i is the same as that found for 
edge-on galaxies in Paper I. For comparison, we show in Fig.2.b the molecular fraction 
for edge-on galaxies taken from Paper I. It is clear that the behavior of f mo \ is the same in 
both figures, regardless of face-on or edge-on. This implies that the deconvolving method 
which we used in Paper I is useful. Fig.2.a and £g.2.b suggest that the phase transition of 
HI and H 2 in galaxies occurs within a narrow range of radius in any galaxies of Sb and Sc 
types. The molecular front is one of the fundamental features in late type disk galaxies. 


3 Gas Phase Transition Theory 

The HI-H 2 gas phase transition theory in galaxies has been developed by Elmegreen (1993, 
hereafter E93), and we apply this theory to the galaxies for which we have presented the 
molecular fraction. We briefly summarize Elmegreen’s phase transition theory. 

H 2 molecules are formed on the surface of dust, while they are destructed by UV 
photons. Crucial parameters for determining / mo \ are P (interstellar pressure), U (UV 
radiation held, noted as j in E93) and Z (metallicity). Here the amount of dust is assumed 
to be proportional to the metallicity Z , and the interstellar pressure P determines the 
density of gas in clouds. It is clear that higher Z and lower U make / mo i higher, and since 
higher density enhances the probability of encounter between HI atoms and dusts, higher 
P also makes / mo i higher. 

Consider a spherical cloud of mass M. Clouds are assumed to be divided into two types 
according to their internal structure; diffuse clouds and self-gravitating clouds. Diffuse 
clouds are assumed to be ideal gas without any internal structure and confined by the 
interstellar pressure. Self-gravitating clouds are assumed to be isothermal gas of spherical 
shape and also assumed to be bounded by interstellar pressure. The critical mass between 
two types is defined by M crit which is the mass at which the average density of diffuse 
clouds becomes equal to the average density of self-gravitating clouds. The critical mass 
is about 1O 4 M 0 in the solar neighborhood, and it of course varies with the interstellar 
pressure P. The temperature of clouds are assumed to be constant through the galactic 
disk, 200K for diffuse clouds and 10K for self-gravitating clouds. These are standard 
values for the clouds in the Milky Way Galaxy. 

Smaller clouds contain no molecules, because UV photons from outside can reach even 
to the center of the cloud and destroy all molecules. On the other hand, larger clouds 
can contain more amount of molecules in the inner part where no UV photon can reach. 
Thus the threshold mass for clouds which can contain molecules changes with P, U and 
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Z. As for clouds which are large enough to contain molecules, we can divide the clond 
into two parts at a certain radius; the inner molecular core and the outer atomic layer. 
The radius at which the cloud is divided into two parts is determined by P, U, and Z, 
as well as by the cloud’s mass. In this way we can calculate the molecular mass of each 
cloud for a given set of P, U, Z and M. 

To obtain the molecular fraction in galactic scale, we have to take into account all 
clouds with various masses. The mass distribution of clouds can be approximated by a 
power law and we take the mass function with the power law index of -1.5 (Combes 1991) 
with the lower cutoff mass of Ml = 1M 0 and the upper cutoff mass of Mu = 1O 6 M 0 , 
following E93. The number of clouds in the mass interval (M, M + dM) is then given by, 

(j) c dM = (f) c0 M~ 15 dM (2) 

Here 0 c o is a constant, and it is determined so that the integration of this mass function 
makes total gas density as below, 

rM u 

p= / </> c0 MM~ 15 dM (3) 

J M l 

Using this mass function we can calculate f mo i in galactic scale by integrating all clouds 
as, 

, Jml (ficoM 1 ' 5 M mo \dM 

/mQl SffifaM-^MdM ’ W 

where M and M mo i are the mass of the cloud and the mass of molecular core of the cloud 
respectively, and Ml and Mu are taken to be 1M 0 and 1O 6 M 0 , as mentioned above. 

Since the resulting f mo \ changes with P, U, and Z, it changes with the galacto-centric 
radius R, because P, U, and Z also change with the radius in real galaxies. 


4 The Molecular Front in a Model Galaxy 

Using the phase transition theory described above, we calculate / mo i as a funtion of 
galacto-centric radius R. To calculate f moh we have to specify P, U and Z as a function 
of radius. We assume that P, U and Z can be given by exponential functions of a radius. 
Although this assumption seems to be rough, several evidences support it; the surface 
brightness of spiral galaxies in the blue band can be well fitted by an exponential function 
(e.g. van der Kruit & Searle 1982; Elmegreen & Elmegreen 1984). Since the distribution 
of the blue light is thought to come from young stars and UV photons come mainly from 
young stars, U can be approximated by an exponential function of the scale length as 
that of the blue band intensity distribution. The metallicity Z is known to decrease 
exponentially with a radius from the observations of many HII regions in our Galaxy and 
nearby face-on spiral galaxies (e.g. Diaz 1989; Belly & Roy 1992; Vila-Costas & Edmunds 
1992). The interstellar pressure P has not been investigated well, but if we assume the 
equipartition between the interstellar hot gas and the random motion of clouds, we can 
estimate the interstellar pressure as follows, 

P ~ Pgas* 7 j (5) 

where p gas is local gas density and a is the velocity of random motion of clouds. Since 
a is about 10 km/s and does not changes so much in our Galaxy, we may assume that 
<j is constant in our Galaxy, and is also so for other galaxies. Then we can obtain the 
interstellar pressure P easily as follows, 


P OC Pgas 


( 6 ) 
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According to Young & Scoville (1982), the H 2 gas density for Sc galaxies is fitted well by 
an exponential function. Since the H 2 gas is dominant in the inner region for which we 
want to explain the behavior of / mol , the total gas density is roughly approximated by an 
exponential function. Then, we may assume that the interstellar pressure P also obeys 
an exponential law. 

Following Elmegreen(1993), we scale P,V and Z with the values of solar neighborhood; 
Pq, Uq, and Z®, and express the three parameters as a function of galacto-centric radius 
as follows, 

(P/Pe) = §-e*P(-R/Rp) (V 

P© 

(U/U e ) = ^eM-R/Ru) (8) 

U Q 

(Z/Z e ) = ^-exp(-R/R z ) (9) 

Z.0 


where Rp , Rp and Rz are effective radii of each parameter, and subscript 0 means the 
value at the galactic center. It is difficult to determine P 0 , Uq, and Z 0 directly, so we 
define Rp s , Rjj s and Rz s , at which (P/P 0 ), (U/U Q ) or (Z/Z®) equals to unity. 

Substituting P(R Ps ) = P 0 , U(Rp s ) = U®, and Z(R Zs ) = Z Q into above equations, P, 
U and Z are written as follows; 

(P/P 0 ) = exp(—(P - R Ps )/Rp), (10) 

(P/P®) = exp(—(P - R Us )/R v ), (11) 

(Z/Z®) = exp(-(P - R Zs )/R z ). (12) 


So, what we have to know are six parameters, Rp and Rp s for calculating P, Rp and Rp s 
for U, and R z and R Zs for Z. 

In real galaxies, necessary parameters for calculating f mo \ must be given independently 
of each other. However, they may be related to each other through the interstellar physics 
and evolution of the galaxy, and it is likely that the scale length of each parameters are 
similar; for example, the optical scale length and the scale length of molecular gas density 
are similar to each other (Young & Scoville 1982), and the optical scale length and the 
scale length of metallicity gradient are similar (Diaz 1989; Zaritsky et al. 1994). In 
order to know the behavior of f mo \ in a general case rather than a certain galaxy, we first 
consider the simplest case; we assume that Rp = Rp = R z and Rp s = Rp s = R Zs , and 
define R eS and R s as follows; 


Res = Rp — Ru — Rz, (13) 

R s = Rp s = Rus = Rzs- (14) 

The former equation means that the scale radii of P, U and Z are identical to each other, 
and the latter means that / mo i = / mo io a t R = R s - 

Fig. 3.a - 3.b 

Under these assumptions, we have calculated the variation of / mol in a model galaxy 
which has an exponential disk. We show the result in Fig.3.a, where we have plotted four 
lines which correspond to the case of R s /R c g =1, 2, 3 and 4. In all cases, the observed 
characteristics of / mo i are reproduced well; / mo i at the galactic center is almost unity, 
/moi at the outer region is almost 0, and the transition between HI and H 2 occurs rather 
abruptly within a narrow range of radius, which we have called the molecular front. It 
is remarkable that the smooth variations of P, U and Z, which are assumed to obey an 
exponential function with the same scale length, indeed reproduce such a drastic change 
°f /moi- We note that the molecular front always appears regardless of the ratio of R s to 
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R e g. This result implies that the molecular front may be one of the fundamental features 
of galaxies which have exponential disks. 

To investigate how each parameter P, U , and Z contributes to / mo i, we consider the 
extreme cases in which one of the three parameters (P, U, Z ) is assumed to be constant 
over the entire disk (equal to the value in the vicinity of the Sun). We show the result 
in Fig.3.b for the case of R s /R e g = 3.5 but one of (P, U, Z) is fixed to be equal to the 
solar value. It is clear from Fig.3.b that f mo \ changes its shape if one of (P, U, Z) is 
assumed to be constant, and that all of the three parameters contribute to f mo \ to some 
extent. However, the most striking feature in Fig.3.b is the case in which Z is assumed to 
be constant. In such a case f mo \ does not reach unity even at the galactic center and the 
slope of / mo i is extremely flat through the entire disk, showing no front-like features at all. 
If we assume that the metallicity Z is constant of 2 Z & , then / mo i increases systematically 
at any radius while keeping the slope nearly identical. On the other hand, there still 
remains front-like feature in the case of constant P or U. This result implies that / mo i 
depends much more on Z than P or U, and that the metallicity gradient is crucial to 
the formation of the molecular front. Since almost all disk galaxies are known to have 
metallicity gradient (Belly & Roy 1992; Vila-Costas & Edmunds 1992; Zaritsky et al. 
1994), we may say that the molecular front is one of the characteristic features resulting 
from the chemical evolution of disk galaxies. We note that the strong dependence of / mo i 
on Z is dne to the facts that the metal absorbs the UV photons and shields H 2 molecules 
from UV destruction as well as it works as a catalysis for the formation of H 2 molecules. 


5 Application to Real Galaxies 

5.1 Face-on Galaxies 

Now we try to reproduce the variation of f mo 1 in real galaxies. We give the distribution of 
P, U, and Z independently and see whether the model with a more realistic distribution of 
P, U and Z affects the behavior of / mo i- This will be done more easily for face-on galaxies 
because the metallicity distributions, which are crucial to reproduce / mo 1 , are available. 
So, we try first to reproduce / mo 1 in face-on galaxies for which we have presented the gas 
distribution and the molecular fraction in section 2. 

Since we have the data of gas density distribution, we can better estimate the interstel¬ 
lar pressure P using the relation (P/P 0 ) ~ (p ga s/Pgaso) rather than using an approximated 
distribution in the form of an exponential function. Furthermore, the surface brightness 
distribution can be obtained from optical observation, and the metallicity and its gra¬ 
dient also can be obtained from optical spectroscopy. Thus, we can give all parameters 
independently and investigate whether the model can describe the actual behavior of f mo 1 
in galaxies. The parameters used here and the references from which the data have been 
taken are summarized in Table 2. 

Fig. 4.a - 4.d 

We show the results from the model calculation in Fig.4.a-4.d, superposed on the 
observed f mo \. As seen in Fig.4, / mol predicted by the model can reproduce the observed 
one. In the case for M101, NGC6946 and IC342, the model predict the molecular front 
clearly, and its shape and location agrees well with the observation although there exists 
slight deviation. In the case for M51, although the observed f mo \ is almost unity within the 
region where the data are available and show no front-like feature, the model reproduce 
the observed f mo 1 well again. The model predicts rather flat f mo 1 , almost constant unity 
in this region and no molecular front at least in this region. This result comes from the 
fact that the metallicity in this region is very high (about 3.5Z Q according to Vila-Costas 

6 Edmunds 1992) and constant up to 6 kpc. As described in section 4, the metallicity 
gradient is crucial to the existence of the molecular front, and we have shown that the 
molecular front would not appear in a galaxy without metallicity gradient (Fig.3.b). The 
case of M51 just corresponds to this case but higher metallicity. It is known that the 
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metallicity decreases exponentially in the outer region of M51 where we do no concern 
here (Vila-Costas & Edmunds 1992). So this galaxy would have the molecular front at 
large radius, (see Kuno et al 1995 for a more detailed study of f mo \ in M51) 

Although the model reproduces the qualitative behavior of f mo \, there still exists resid¬ 
ual to be explained. We think it comes from both the observational error and the model. 
As mentioned above, we give U and Z by exponential functions, but it is likely that the 
true distribution of U or Z deviates from an exponential function (locally and/or glob¬ 
ally.) For instance, the metallicity gradient in the inner region of M51 and M101 differs 
from that in the outer region (Vila-Costas and Edmunds 1992). This result implies that 
the metallicity gradient in a form of exponential function is not universal, but sometimes 
not appropriate. So a more precise observation of the distribution of the metallicity is 
necessary for further argument. In addition to this, the quantitative disagreement of / mo i 
also comes from the use of the optical scale length at the blue band instead of the U 
band. In blue band, the contribution of evolved giants to the luminosity profile is not 
entirely negligible, and the surface brightness may not trace correctly the distribution of 
young OB starts, while this is not the case in U-band in which young OB stars dominate 
the luminosity. So we need the data of U band surface brightness distribution with high 
quality for further improvement. 

On the other hand, we also think that the model itself contains some uncertainties. 
The phase transition theory by Elmegreen (1993) includes several parameters which have 
not been determined precisely yet; for example, the power low index of the clouds’ mass 
function, lower and upper limit masses for the clouds’ mass function, the temperature of 
clouds, and so on. Their values used in the theory are thought to be common now, but 
they might be changed by further observation. 

5.2 Edge-on Galaxies 

For edge-on galaxies, for which f mo \ have been obtained in Paper I, there are some dif¬ 
ficulties in obtaining the accurate values of necessary parameters, especially parameters 
which are related to the metallicity (the metallicity and its gradient) because they have 
been measured only in nearby face-on spiral galaxies. Here we estimate the metallic¬ 
ity by using empirical relations, since the metallicity gradient normalized by the optical 
scale length of the disk is almost constant through galaxies (Diaz 1992; Zaristky et al. 
1994), and the characteristic abundance of galaxies log(O/H) 08i?s , which is defined as 
the metallicity at the radius of 0.8xRb ( Rb ; the optical scale length at the B band) is 
correlated to the mass of galaxies (Zaritsky et al. 1994). These empirical relations are 
shown in Fig.5.a and Fig.5.b, in which the data are taken from Zaritsky et al. (1994). 

Fig.5.a - 5.b 

We have obtained the empirical relation between the characteristic abundance and the 
rotation velocity of galaxies, and the relation between the metallicity gradient and the 
rotation velocity by applying a least square fit. The resulting relations are expressed as 
follows; 

[12 + log(O/H )] 0 _ 8Rb = 8.53 ± 0.17 + 0.0244V rot , (15) 

A[log(0/H)]/A B = -0.196 ± 0.105 + 6.81 x 10~'V rot . (16) 

The first relation shows that more massive galaxy has larger characteristic abundance. 
This tendency is consistent with the results obtained by Roberts & Haynes (1994). In 
the second relation, the dependence of the gradient on the rotation velocity is small and 
the gradient is almost constant. The constant value of -0.196 is consistent with the result 
by Diaz (1989) when we consider the scatter. Since these relations are not unique and 
have some scatter, we cannot give a unique value of Rz or Rz s from these empirical 
relation. We fix the metallicity gradient to be equal to the moderate value in the second 
relation, while two possible characteristic abundances are considered; the maximum and 



the minimum characteristic abundance which correspond to ±lcr in the first relation. 
So the metallicity gradients in two cases are the same but the total amount of metal is 
different from each other. 

The data of surface brightness distribution for NGC 3079 have been taken from Shaw 
et al. (1993), and the data for other galaxies from van der Kruit & Searle (1982). The 
interstellar pressure P is estimated from the gas distribution, in the same way done in 
the former subsection. The parameters for these galaxies are summarized in Table 3. 

Fig. 6.a - 6.c 

We show thus obtained / mo i for NGC 891, NGC 3079, and NGC 4565 in Fig.6.a-6.c, 
superposed on the observed / mo i which have been obtained in Paper I. We exclude NGC 
5907 because the HI data for this galaxy has not been obtained over the hole disk and it 
might cause for substaintail error in estimating the interstellar pressure P, especially in 
the outer galaxy where HI gas is dominant. In Fig.6.a, we show the results for NGC 891. 
The model reproduces the variation of / mo i well again, as seen for the face-on galaxies. 
The observed / mo i lies between two plots and it is slightly nearer to the case of the higher 
characteristic abundance. So the characteristic abundance in this galaxy may be really 
higher than the moderate value in the relation obtained above. For the case of NGC 3079, 
the observed / mo i is also in the possible range predicted by the model and the model can 
reproduce f mo \ well. The scatter of the empirical relation causes less difference of f mo \ than 
that for NGC 891, because the molecular gas is highly concentrated in the central region 
of this galaxy (Paper I; Sofue & Irwin 1992) and the effect of the interstellar pressure P 
must be larger than the effect of the scatter of the Z. It is remarkable that the present 
model is applicable to such a peculiar galaxy with nuclear activity. This implies that the 
interstellar physics in the disk of active galaxy such as NGC 3079 must be the same as 
the normal disk galaxies. In Fig.6.c we show the models for NGC 4565. The model can 
explain the qualitative characteristic of f mo i for NGC 4565, including the sudden decrease 
at R~3 kpc and the enhancement at R ~ 5kpc. The f mo \ calculated by the model for 
NGC 4565 is smaller than the observed / mo i throughout the disk, but if we give a slightly 
steeper gradient, it is possible to reproduce the observed / mo i, because f mo \ is so sensitive 
to Z , and its value may lie in the possible range in Fig.5.b. 

The observed f mo \ for edge-on galaxies has been reproduced well using moderate metal¬ 
licity gradient which is determined by the empirical relation. In the former subsection we 
have already seen that the observed / mo i for face-on galaxies has also been explained well 
using the set of (P, U, Z) determined observationally. From these results we stress that 
the model described in the present paper works well to explain the phase transition of HI 
and H 2 gases in galaxies. 


6 On the Relation to the Chemical Evolution of the Galaxy 

As described in Section 3.2, f mo \ depends strongly on the metallicity Z and the metallicity 
gradient is crucial to the formation of the molecular front. These results imply that the gas 
distribution observed at present is tightly related to the history of the chemical evolution 
of a galaxy. Since the disk of spiral galaxy was virtually metal-free at the earliest stage 
of the evolution and the metallicity has been increasing monotonically with time, / mo i 
at any place in the disk must have changed with time. Then, the question is how f mo \ 
has varied with time in a galaxy, and when and how the molecular front appeared and 
evolved since then. Since there was no metallicity gradient in the disk at the early stage, 
the molecular front did not exist at that time. It must have formed with the growth of 
the metallicity gradient, but when and how it has been formed is still a matter of debate 
because it is not known yet how the metallicity gradient formed in a disk galaxy. 

Once the molecular front is formed, it would not disappear unless the metallicity 
gradient would vanish. However, the molecular front would not stay at the same position, 
instead it should change its shape and location as the chemical evolution proceeds in a 
galaxy. For example, suppose that the metallicity increases in global while keeping the 
slope of the gradient constant. In this case, the molecular front advances outwards with 
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time without changing its shape; this is seen in Fig.3.a. The most chemically-evolved case 
corresponds to the line with R s /R e s=4, while the least enriched case corresponds to the 
case of R s /R e ff = 1. In real galaxies the slope of the metallicity gradient may also change. 
However, there is no doubt that the metallicity has increased monotonically, so we think 
the molecular front may advance outwards with time. 

Though it is still controversial whether the star formation rate is more tightly related to 
HI gas or to H 2 gas, it is interesting to consider the effect of / mol on the evolution of galaxies 
through star formation. We here assume that the star formation rate is determined by 
the amount of H 2 gas, not by HI gas. More precisely, we should assume that the star 
formation rate is determined by the amount of molecular gravitating clouds rather than 
total amount of molecular clouds which include both gravitating and diffuse molecular 
clouds, because star formation would rarely occur in diffuse clouds. However, as far as our 
model is concerned, more than 80% of H 2 gas is made up from self-gravitating clouds in 
any case, so it is a good approximation to assume that star formation rate is determined 
by the amount of H 2 . If we assume a Schmidt’s law, SFR oc n(H 2 ) Q , then f mo \ represents 
somehow the star formation efficiency from total gas. So the variation of / mo \ through the 
disk manifests the star formation efficiency in the disk; star formation efficiency is high 
in the inner disk and many stars are forming or will form soon there, while few stars are 
forming in the outer region. So, the molecular front divides the disk into two parts; the 
inner star forming disk and the outer disk where stars are rarely formed. This is indeed 
observed in NGC 891 (see Sofue & Nakai 1993). 

Once stars are formed under a condition of certain / mo i, newly formed OB stars radiate 
UV photons and soon explode as supernovae, ejecting heavy elements. So some 10 6 ~ 7 
years after the star formation, the physical condition of ISM must be changed by the 
newly formed stars. This implies the change of (P,U,Z), and so the change of / mol and 
star formation efficiency. A galaxy must evolve with such a cycle; a given set of (P,U,Z) 
determine f mo \ and so the star formation efficiency, and then newly formed stars under 
such condition of star formation changes the physical condition; P , U and Z , leading to 
different f mo \ and star formation efficiency. 

From such a view of galaxy evolution, it is interesting to consider the star formation 
history at the earliest stage of evolution of a galaxy. The first generation of stars must 
have been formed from the unenriched molecular cloud in the central region of a galaxy, 
where H 2 molecules had been formed by a direct triple encounter of H atoms. (We do 
not take into account this process for calculation of / mo i because the formation of H 2 by 
dust is much more efficient than direct encounter, but in case of no metal, H 2 molecules 
must have been formed by such process.) Once stars formed, they emit UV radiation and 
destruct molecules quite efficiently compared with the molecular clouds at the present 
epoch because of the lack of absorbing dust. So star formation efficiency soon decreases 
with the star formation, and so only few stars are formed after the first star formation. 
Then some 10 6 ~' years later all of OB stars of the first generation explode as supernovae, 
ejecting heavy elements, and this makes U lower and Z higher, and so f mo \ higher, leading 
to the formation of second generation of stars. A galaxy at very early stage must have 
evolved with such a drastic change of / mol until the gas has been enriched enough to 
protect clouds from the UV radiation. 

In summary, in order to investigate the evolution of galaxy more precisely, the atomic- 
to-molecular phase transition process is important through the star formation. So, it 
would be of vital importance to construct models which include such a process instead 
of the simple star formation low such a Schmidt law, in which the star formation rate is 
simply assumed to proportional to the total amount of gas. 
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Figure captions 

Fig.l.a The surface mass distribution of Hf and H 2 gases in M51. The amount of H 2 
gas is calculated from the intensity of CO line, assuming the constant conversion factor 
of 3xl0 20 . The references from which the data have been taken are listed in Table 1. 

Fig.l.b Same as Fig.l.a, but for M101. 

Fig.l.c Same as Fig.l.a, but for NGC 6946. 

Fig.l.d Same as Fig.l.a, but for 1C 342. 

Fig.2.a The molecular fraction / mo i in face-on galaxies. 

Fig.2.b The molecular fraction f mo \ for edge-on galaxies, which has been obtained in 
Paper I. 

Fig.3.a The molecular fraction f mo \ for a model galaxies in which all of P, U. and Z 
obey an exponential function with same scale length of R c g. These four lines correspond 
to the case that the ratio of equals to 1, 2, 3 and 4. In all cases / mo i changes 

rather drastically at a certain region, which we call molecular front. 

Fig.3.b Models without gradient of one of three parameters P, U, and Z. The solid 
line corresponds to the model same as in Fig.l with R/R c r = 3.5, and other three lines 
correspond to the cases of constant P, U, or Z. It is remarkable that the metallicity Z is 
the most dominant parameter. 

Fig.4.a f mo i for M51 by the model, superposed on the observed f mo \. The necessary 
parameters are given observationally and are listed in Table 2. The vertical bar indicates 
typical observational error of / mo i, and the error in the inner region is less than that. 

Fig.4.b same as Fig.4.a, but for M101. 

Fig.4.c same as Fig.4.a, but for NGC 6946. 

Fig.4.d same as Fig.4.a, but for IC 342 

Fig.5.a The empirical relation between the characteristic abundance and the rotation 
velocity of galaxies. Characteristic abundance is defined as the abundance at R=0.8i?^ 
where R# is the scale length of the optical disk. The solid line shows the empirical relation 
given by least-square fit, and two dashed lines show the scatter of la. 

Fig.5.b The empirical relation between the metallicity gradient normalized by the 
scale length of the optical disk and the rotation velocity. Note that the gradient is almost 
constant through galaxies. The notation of solid and dashed lines are same as Fig.5.a. 

Fig. 6 . a f mo i calculated for NGC 891, superposed on observed / mo i from Paper I. The 
dotted line corresponds to the case with the moderate Rz and Rz s , and the heavy solid 
line corresponds to the case of the maximum R z and Rz s (la). The observed f mo \ (thin 
line) lies between two plots, but is much nearer to the case of of the maximum Rz and 
Rzs 

Fig.6.b Same as Fig.6.a, but for NGC 3079. 

Fig.6.c Same as Fig.6.a, but for NGC 4565. 
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Table 1 : Assumed distances and references for four face-on galaxies. 


galaxy 

Distance 

references for distance 

references for HI and CO 

M 51 

9.6 Mpc 

a 

c,d 

M 101 

4.1 Mpc 

b 

e 

NGC 6946 

5.9 Mpc 

b 

f 

IC 342 

2.7 Mpc 

b 

f 


a. Sandage & Tammann (1974). 

b. Vila-Costas & Edmunds (1992). 

c. Rots et al. (1990). 

d. Nakai et al. (1994). 

e. Kenney et al. (1991). 

f. Young & Scovillc (1982) 
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Table 2 : The parameters for four face-on galaxies. 


galaxy 

E> a 

Ku 

D a 

ttUs 

Rz 

Rzs 

M 51 

5.2 

11.5 

— 

_b 

M 101 

5.2 

12.3 

3.3 

5.2 c 

NGC 6946 

4.0 

6.7 

4.9 

6.5 d 

IC 342 

5.5 

2.9 

5.0 

5.9 c 


a. Elmegreen & Elmegreen (1984). 

b. The metallicity in the inner region of M51 is constant, equal to 3.5 Z Q (Vila-Costas 
& Edmunds 1992). 

c. Vila-Costas & Edmunds (1992). 

d. Belly & Roy (1992). 
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Table 3 : The parameters for three edge-on galaxy. 


galaxy 

Distance a 

Ru 

Rus 

NGC 891 

8.9 Mpc 

4.6 kpc 

5.0 kpc b 

NGC 3079 

15.0 Mpc 

4.7 kpc 

7.3 kpc c 

NGC 4565 

10.2 Mpc 

5.5 kpc 

9.5 kpc b 


a. Paper I. 

b. van der Krnit & Searle (1982). 

c. Shaw et al. (1993). 
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